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Neutrino telescopes may have the potential to detect the quasi-stable staus predicted by some 
supersymmetric models. Detection depends on stau electromagnetic energy loss and weak interac- 
tions. We present results for the weak interactions contribution to the energy loss of high energy 
staus as they pass through rock. We show that the neutral current weak interaction contribution is 
much smaller than photonuclear energy loss, however, the charged current contribution may become 
dominant process above an energy of ~ 10 9 GeV. As a consequence, the stau range may be reduced 
above ~ 10 9 GeV as compared to the range neglecting weak interactions. We contrast this with the 
case of tau range, which is barely changed with the inclusion of charged current interactions. 



I. INTRODUCTION 

Interactions of very high energy neutrinos with nucle- 
ons as they traverse the Earth are ideal probes of physics 
beyond the Standard Model High energy neutrinos 
originate in interactions of high energy cosmic rays with 
microwave background photons (cosmogenic neutrinos), 
or they might be produced in astrophysical sources such 
as active galactic nuclei and gamma ray bursts Q • These 
neutrinos do not interact on their way to the Earth, and 
they arrive undeflected by magnetic fields. Once they 
reach the Earth, they interact with nucleons in the Earth, 
or possibly in the atmosphere. In some supersymmet- 
ric models, neutrino interactions in Earth could produce 
heavy supersymmetric particles that decay into quasi- 
stable sleptons. 

In the weak scale supersymmetric models with the su- 
persymmetry breaking scale larger than 5 x 10 6 GeV, 
the next-to-lightest particle (NLP) is a charged slepton 
(stau) which eventually decays into the stable lightest su- 
persymmetric particle (LSP), the gravitino [||. Due to 
very long lifetimes, staus may travel thousands of kilome- 
ters through the Earth without decaying. Recently it was 
proposed that a direct way of probing the SUSY breaking 
scale in weak scale supersymmetry models would be to 
detect pairs of charged tracks in neutrino detectors, such 
as IceCube, from staus resulting from neutrino-nucleon 
interactions producing heavier supersymmetric particles 
[Hi- This has been further explored in Refs. [j| and [6[. 

Because of the small cross section, the production of 
staus from downward neutrinos is negligible in compari- 
son with the background from the standard model pro- 
cesses. However, upward neutrinos producing staus could 
potentially be detectable because of the effective detector 
volume that is enlarged by the long range of the stau. 

Detection of staus produced in neutrino-nucleon inter- 
actions in Earth depends strongly on the stau lifetime 
and range. Thus it is crucial to determine the energy 
loss and effective range of the high energy stau as it tra- 
verses the Earth. The average electromagnetic energy 



loss of a particle which traverses a distance X is given by 



dE 



(1) 



where E is the lepton energy, a represents the ionization 
energy loss, and (3 is the radiative energy loss. To first 
approximation, p scales inversely with stau mass. 

In Ref. 0, we evaluated the electromagnetic energy 
loss of scalar leptons more quantitatively. We showed 
that the photonuclear interaction gives the largest con- 
tribution to (3 for stau energies between 10 6 — 10 12 GeV. 
The range determined by electromagnetic interactions is 
of order 10 4 km.w.e. for stau masses of a few hundred 
GeV. Interaction lengths from charged current weak in- 
teractions are of the same order of magnitude. In this 
paper, we evaluate the energy and mass dependence of 
the stau weak interaction energy loss and attenuation, 
and we show its relevance to the stau range. 

Weak interaction cross sections, because of the mas- 
sive vector boson propagators, have a different mass de- 
pendence than the electromagnetic energy loss parameter 
(3. Depending on the details of the supersymmetric cou- 
plings, we find that weak effects may dominate the stau 
range. 

In Sec. II, we present results for neutral and charged 
current cross sections and energy loss for staus and dis- 
cuss its mass dependence. We make comparisons with 
the lepton case. We show the tau and stau range includ- 
ing weak interaction processes using a one-dimensional 
Monte Carlo evaluation described in Ref. Q in Sec. III. 
The weak interaction effects in r propagation through the 
Earth are determined to be negligible for the energies 
considered here - up to 10 12 GeV. Our conclusions for 
the tau range differ from estimates based on characteris- 
tic distance scales for the tau @ . We discuss implications 
of the stau range including maximal charged current in- 
teractions for IceCube and higher energy measurements 
such as by ANITA in Sec. IV. 
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II. STAU ENERGY LOSS 

The energy loss parameter (3 has contributions from a 
variety of processes: 

where y is the fraction of lepton energy loss in the radia- 
tive interaction, 



E-E' 



(3) 



for final stau energy E' . The superscript i denotes 
bremsstrahlung (brem) ficl ITU , pair production (pair) 
|T^ |. photonuclear (nuc) |8l.ll3| and weak (NC) processes 
for interactions of the initial particle with a target nu- 
cleus. Avogadro's number is Na and the atomic mass 
number of the target nucleus is A. 

At low energies, where j3E <C a, either the lifetime or 
ionization energy loss determines the stau range, which 
scales linearly with energy 0] . The ionization energy loss 
parameter a is nearl y co nstant as a function of mass of 
the particle, namely [lj] 



2 x l(T d GeVcmVg 



(4) 



At energies above 10 6 GeV, energy loss of leptons and 
staus is dominated by the electromagnetic radiative pro- 
cesses. For staus, weak interactions may also become 
important. We review next the weak interaction cross 
sections for scalars. 



A. Weak Interaction Cross Sections 

Neutral current and charged current cross sections are 
relevant in two different ways. Neutral current interac- 
tions only shift the stau energy. The neutral current 
interactions can be incorporated into (3. By contrast, the 
charged current interactions remove the stau, in the pro- 
cess producing a sneutrino. The sneutrino then decays, 
presumably to another stau. We do not include stau 
regeneration because of the several steps decreasing en- 
ergy. Neutral current and charged current cross sections 
of staus are shown below, as well as the results for the 
cross sections for taus. 



1. Neutral Current Cross Sections 

The neutral current cross section that describes the 
interactions of charged sleptons (f ) with nucleons via ex- 
change of Z° boson is given by 



d 2 a NC {fN) 
dxdy 
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(5) 



The parton fractional momentum is x and y is the frac- 
tion of slepton energy loss. The quantities a/ and (3f are 
the couplings of staus to gauge bosons [l5j] 



OLf = — ^3 tan 9w — cot 



Ow 



and 



/?/ = -( tan 9 W + cot Ow 



(6) 



(7) 



The scalar partner of the right-handed tau may not be 
a mass eigenstate. The angle Of parameterizes the mix- 
ing between scalar partners of the right-handed and left- 
handed tau, where sin Of = means that the mass eigen- 
state quasi-stable stau is purely made of the partner of 
the right-handed tau. In principle, smOf need not equal 
zero. We take smOf — 1 for the neutral-current process 
in the figures below to evaluate the maximal effect in 
the charged current case, since W's couple only to left- 
handed fermions and their scalar partners. The range of 
(ctf + /3fCos20f) 2 is such that 

< Uaf+(3fcos20f) 2 /(af+(3f) 2 =r NC ] < 1.38 . (8) 

For the neutral current structure functions, we have 
taken 2xFi — F2 here and for charged current interac- 
tions. For neutral currents, 



Fx C = ^i+a 2 )[ qi (x,Q 2 )+^,Q 2 )] (9) 



1 



with 



Vi = T 3 - 2ei sin 2 W 



(10) 
(11) 



for weak isospin assignments T3 = ±1/2 and electric 
charge e^. We use CTEQ6 parton distribution functions 
fl6l ] with a power law extrapolation of these distributions 
for x < 10~ 6 of the form x~ Xi , where i denotes quark or 
antiquark flavor [Tt], [ill]. The values we use for \ are 
given by: 

uu, dd ss, bb cc 

A, -0.0276 -lnQ + 0.1784 \ u + 0.0054 A u + 0.0094. 

The kinematic limits on the variables of integration, y 
and Q 2 (for small y) are given by 

^ < Q 2 < 4E 2 (l~y)- m ^~y? 



f(M + m n ) 2 - M 2 ) 
2ME 



< y< 1 



nif 



(12) 
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In Fig. 1 we show the neutral-current cross sections 
for stau masses of 50 GeV and 250 GeV, for sindf = 1. 
We also show the muon, tau and neutrino neutral current 
cross sections for comparison. We note that the cross sec- 
tion for taus is almost indistiguishablc from the muons 
because the masses of taus and muons are small com- 
pared to the energy considered. The stau cross sections 
have weak mf dependence. 

We see that the stau NC cross section is almost an 
order of magnitude smaller than the neutrino case for 
sm0 f = 1 and m? = 250 GeV at E = 10 6 GeV. The 
difference comes from the couplings as well as the y de- 
pendence of the differential cross section. The ratio of 
neutrino to stau neutral-current couplings, including spin 
averaging, yields about a factor of 2. At small y, the 
differential cross section, da/dy, for neutrinos is larger 
by about a factor of 5 than for the stau. The aver- 
age y for staus decreases with increasing stau mass from 
< y >« 0.13 for m f = 50 GeV to < y >« 0.06 for m? = 
250 Ge 



and we use the kinematic limits 
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FIG. 1: Neutral-current cross sections for neutrino, tau, 
muon, and stau for mf = 50 GeV and rtif — 250 GeV, with 
sin (9/ = 1. 



2. Charged Current Cross Sections 



m%y 
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(15) 



where mf represents the mass of the incoming stau and 
mo represents the outgoing sneutrino. For the stau pro- 
cess we take — m? = 50 GeV, with m? as 50 GeV and 
250 GeV. 

In Fig. 2 we show the charged-current (CC) cross sec- 
tions for the stau with mass 50 GeV and 250 GeV, and 
for the tau, muon and neutrino for comparison. We note 
again that the cross section for stau has a weak mass 
dependence. The cross sections for taus and muons are 
indistinguishable due to the small masses relative to the 
energies we consider. The charged lepton CC cross sec- 
tion is about a factor of 2 smaller than the neutrino case, 
due to the spin averaging. The energy dependence of the 
stau cross section is stronger than for the tau and muon. 




— X. U- 

- . m st = 50 GeV 

. - m = 250 GeV 



10 10 



10° 10' 10'" 10" 

E [GeV] 



The charged-current cross section of staus with nucle- 
ons via exchange of charged boson is given by 



d 2 a cc (fN) 
dxdy 



G 2 F 



( \ 



y 



2ME 



sin 2 Of ME 
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(13) 



Since we are interested in considering the upper limit on 
the stau cross section, we take sin 2 9 f = 1, however, the 
value of s'm9f is unknown. For the charged current, the 
conventional normalization of the structure functions is 



^3 



2fe(z,Q 2 



q 3 (x,Q 2 )} 

- qj(x,Q 2 )} 



(14) 



FIG. 2: Charged-current cross sections for neutrino, tau, 
muon, and stau for mf = 50 GeV and mf = 250 GeV, with 
sinOf = 1. 



B. Application to Energy Loss 

The energy loss (3 for neutral current interactions can 
be found at a fixed initial charged lepton or stau energy 
from Eq. ([2]). In Fig. 3 we show our results for (3 NC for 
muon, tau and stau with masses 50 GeV and 250 GeV, 
with sin Of = 1. The values of fl NC for the muon and tau 
are very close in value, as they were for the cross sections. 
Using the relation 



p ~ N(y)a(E) , 
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FIG. 3: Neutral-current f3 for tau, muon, and stau for mf — 
50 GeV and m f = 250 GeV, with sinfl/ = 1. 



current interactions become significant at higher energies, 
in contrast to the neutral current interactions. For a stau 
mass of 250 GeV and sin#/ = 1, this corresponds to 
energies higher than about 4 x 10 9 GeV. For the lighter 
50 GeV stau mass, the charged current process does not 
contribute significantly for energies up to 10 12 GeV. The 
opportunity for charged current interactions to become 
relevant comes from the fact that the charged current 
cross section is less sensitive to the stau mass than the 
photonuclear energy loss parameter f3. 

A comparison of Figs. 5 and 3 verifies our assertion 
that the weak neutral current contribution to (3 for staus 
is not important. This is also true for tau energy loss. 



we estimate the average value for y to be about 0.2 for 
tau and muon neutral-current interaction. We see that 
the stau case is smaller than tau by about a factor of 
3 at 10 6 GeV when we use a stau mass of 50 GeV. We 
show below that the photonuclear contribution to (3 is at 
least an order of magnitude larger than (3 N for staus. 
For muons and taus, the photonuclear (3 nuc is larger than 
(3 NC by a factor - 10 3 - 10 4 . 

In Fig. 4 we show the mass dependence of (3 NC for 
different initial energies. We note that the mass depen- 
dence is weaker than 1 / for ra? < 200 GeV. For masses 
larger than 200 GeV, 1 / rrif scaling works reasonably well. 
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FIG. 4: Mass dependence of neutral-current f3 for stau with 
a stau mass range of rrif — 50 GeV to rrif — 500 GeV, with 
sin (9/ = 1. 

In charged-current processes the initial and final states 
are different so we use Na , the inverse of the effective 
interaction length, instead of (3. We show in Fig. 5 the 
results for Na together with j3 nuc . These figures are 
plotted using sin Of = 1 and for a sneutrino mass 50 
GeV more massive than the stau. For the sneutrino mass 
between 5 GeV and 150 GeV more than the stau mass, 
the cross section changes by a factor of ~ 1/2 — 2. In the 
following, we show only a 50 GeV mass difference. 

Comparing the scales associated with weak interaction 
and electromagnetic processes shows that the charged- 




E [GeV] 



FIG. 5: Charged-current interaction length for tau, muon, 
and stau for m,f = 50 GeV and mf = 250 GeV, with 
sin (9/ = 1. The sneutrino mass is mf + 50 GeV. Also shown 
are photonuclear energy loss parameters for the two stau 
masses. 



III. STAU AND TAU RANGE 

Because (3 NC <C j3 nuc for staus, taus and muons, we ne- 
glect neutral current interactions in our evaluation of the 
range. For staus, the charged current interaction length 
is roughly comparable to 1//3 for some stau masses, so 
we include the charged current interactions in our Monte 
Carlo evaluation of the particle range. Details of the 
Monte Carlo evaluation appear in Refs. [|[ and @- The 
Monte Carlo computer program computes survival prob- 
abilities P(E, Eo, X') for an particle incident with energy 
E which survives a distance X' with E > Eq. The range 
is defined by 

X(E,E ) = J dX'P(E,E ,X') . (16) 

We have taken Eq — 10 3 GeV for the figures shown here 
for staus. 
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FIG. 6: Characteristic distances in kilometers water equiva- 
lent units (dashed) and ranges (solid) for staus in rock, for 
mr = 150 GeV, sin6>/ = and 1, and Vf = 10 7 GeV. The 
minimum stau energy is Eo = 10 3 GeV. The sneutrino mass 
is nif + 50 GeV. 



In Fig. 6, we show characteristic distances associated 
with stau interactions in rock. The curves are evaluated 
for rnf = 150 GeV and the decay parameter F 1 / 2 = 10 7 
GeV, where the lifetime is determined by 



F 



100 GeV 



10 14 GeV' 



10 km 



(17) 



The lifetime is not relevant for this energy range for 
m f = 150 GeV and F 1 ' 2 = 10 7 GeV: EcT P /m f ~ 3 • 10 4 
kmwe for E = 10 6 GeV. A distance which is relevant, 
which also grows with energy, depends on the ionization 
energy loss parameter a through d ~ E/a, shown in 
the figure. We also show the charged current interaction 
length (Na 00 )" 1 and the distance characterized by 
the electromagnetic energy loss parameter. We show the 
range for 150 GeV staus with no charged current contri- 
butions (sin#/ = 0) and with maximal charged current 
contributions (sin (9/ = 1). At low energies, the ioniza- 
tion energy loss dominates, but for E ~ 10 8 GeV, the 
charged current interaction dominates the evaluation of 
the range if smOf = 1. The range does not precisely 
equal the charged current interaction length because the 
electromagnetic energy loss is still a factor, shifting the 
initial stau energy to lower energies. 

Fig. 7 shows the stau ranges in rock for rrif = 150 and 
250 GeV, again for minimum stau energy of 10 3 GeV 
and F 1 / 2 — 10 7 GeV. The upper curves show the range 
when charged current interactions are vanishing, while 
the lower curves have maximal charged current interac- 
tions. 

Finally, in Fig. 8 we show the characteristic distances 
for taus in rock and the tau range. The lifetime governs 
the range at low energies, while electromagnetic energy 
loss dominates at high energies. Because the charged cur- 
rent interaction length is small compared to (/3p) _1 , the 



tau range changes very little with charged current inter- 
actions included. We do not find a decrease in the range 
in rock or water near E = 10 12 GeV as suggested in Ref. 
Q. This is due to the fact that using just the scales, 
e.g., (/3p) _1 or the CC interaction length, is insufficient 
to accurately compute the range. We have directly eval- 
uated the tau range via Eq. (16), where the probability 
includes stochastic effects in the tau propagation. 
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FIG. 7: Range of stau in rock, for m? = 150 and 250 GeV, 
sin (9/ = and 1. The lifetime is here governed by y/~F = 10 7 
GeV, and the minimum stau energy is Eo = 10 3 GeV. The 
sneutrino mass is mf + 50 GeV. 
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FIG. 8: Characteristic distances (dashed) and tau range in 
rock (solid) in km. By including the charged current inter- 
action length, the tau range is unchanged on the scale of the 
figure. 



IV. CONCLUSIONS 

We have shown that weak interactions have the poten- 
tial to play an important role in stau detection by neu- 
trino telescopes, however, the effect is strongly energy 
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dependent. Our results are based on a sneutrino mass 
50 GeV more massive than the stau, and we considered 
sinOf = 1. 

Recent work on stau signals [5j, 16J has focused on the 
IceCube detector. The stau pair event rate depends on 
the effective volume for creation of a pair of staus. The 
effective volume scales approximately with the stau range 
which depends on energy. 

Refs. 0, Q show that the signal events come from 
energies fairly near the threshold for squark production 
because of the falling neutrino fluxes. In these studies, 
an incident E~ 2 neutrino flux is assumed. The threshold 
depends on squark masses: for niq — 300 — 900 GeV, 
the energy threshold is between a fewxlO 5 — 10 6 GeV. 
For E = 10 6 GeV, even a maximal charged current cross 
sections with mo — TUf + 50 GeV doesn't affect the stau 
energy range dramatically. For rrif — 150 GeV, maximal 
CC interactions reduce the range by 22% at 10 6 GeV, 
and by 31% at 10 7 GeV. The effect is less pronounced 
for m? = 250 GeV. When sin 9 f = 1 for m f = 250 GeV 
and F 1 / 2 = 10 7 GeV, the range is reduced by 9% and 
21% for E = 10 6 and 10 7 GeV, respectively. Based on 
these reductions in the range, the thresholds in the 10 5 — 
10 6 GeV energy range and steeply falling fluxes, event 
rate estimates without including weak interactions are 
reasonably reliable. 

Efforts to try to detect staus with higher energy thresh- 
olds are potentially strongly influenced by charged cur- 
rent interactions, where the range can be as much as 
two orders of magnitude shorter than the range evalu- 
ated without charged current interactions. A detector 
such as the Antarctic Impulse Transient Array (ANITA) 
[l9| is sensitive to stau energies larger than ~ 10 s GeV. 
Designed to use a radio antenna suspended by a balloon 
~ 37 km over the south pole ice, the primary goal of 
the experiment is to detect cosmogenic neutrinos inci- 
dent just below the horizon which interact with the ice. 
The goal is to detect the radio Cherenkov signal produced 
by neutrino weak interactions which refracts on its way 
out of the ice. Staus would also make a signal by weakly 
interacting in the ice. 

Weak interactions play a role for ANITA signals in 
two ways. As noted above, the most important feature 
of weak interactions is to produce the signal itself. For 
electromagnetic interactions, only a small amount of en- 
ergy will be deposited in the shower over the area of 
~ 10 6 km 2 . Neutral current interactions have a larger en- 
ergy deposition necessary for detection. Charged current 
interactions have the largest energy deposition. Only 



neutral current processes will contribute to the signal if 
smQf — 0, while charged current processes contribute 
with increasing values of sin 9 f . The maximum con- 
tribution from charged current interactions occur when 
sin 9 f = 1. For this maximum value of the mixing angle, 
10, meaning that the probability for inter- 
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actions of staus in the ice is increased by a factor of 10 
over the case of sin#/ = 0. 

A second effect for ANITA signals is that charged cur- 
rent weak interactions may attenuate the stau flux in 
transit to the ice in view of the detector. By incorpo- 
rating the attenuation of the GZK neutrino flux [gcj as 
it traverses the Earth, the stau production cross section 
4] and the stau interactions on the way to the detector, 
we find that the stau flux attenuation at 10 9 GeV ranges 
from 1 at 0° to ~ 1/3 at a 10° angle measured relative to 



the horizon [2l|. For energies of 10 10 GeV, attenuation 
causes the stau flux to be lowered by a factor of ~ 1/10 
with maximal charged current interactions for an angle 
of 10° below the horizon. The factor of 10 increase in the 
signal due to the stau CC interactions in ice is sufficient 
to compensate for the attenuation as well as enhance the 
signal when the full range of energy is considered. 

To summarize, weak interaction effects are small in the 
range of energies relevant to the IceCube detector. Re- 
cent event rate estimates [3, [H, Q for IceCube will be 
reduced by less than - 30% for E = 10 6 - 10 7 GeV 
by including maximal weak interactions. The potential 
for observing staus at higher energies, for example, by 
the ANITA detector, is enhanced by maximal weak in- 
teractions. The enhancement is due to the high energy 
charged current cross section which becomes increasingly 
important as the stau energy increases above 10 s GeV. 
A detailed investigation of the ANITA signal from stau 
pairs, and for the proposed ARIANNA [22| telescope is 
in progress [2l|. Signals of staus from several energy 
regimes may impose constraints on this class of super- 
symmetry models with quasi-stable staus in the future. 
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